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Abstract 

(I) We exhibit a set of 23 points in the plane that has dilation at least 1.4308, improving the 
previous best lower bound of 1.4161 for the worst-case dilation of plane spanners. 

(II) For every n > 13, there exists an n-element point set S such that the degree 3 dilation 
of S equals 1 + -s/S = 2.7321.. . in the domain of plane geometric spanners. In the same domain, 
we show that for every n > 6, there exists a an n-element point set S such that the degree 4 

dilation of S equals 1 -|- y^( 5 — •\/5)/2 = 2.1755 ... The previous best lower bound of 1.4161 holds 
for any degree. 

(III) For every n > 6, there exists an n-element point set S such that the stretch factor of 
the greedy triangulation of S is at least 2.0268. 

Keywords: geometric graph, plane spanner, vertex dilation, stretch factor. 


1 Introduction 

Given a set of points P in the Euclidean plane, a geometric graph on P is a weighted graph G = {V, E) 
where V = P and an edge uv € E is the line segment with endpoints u,v € V weighted by the 
Euclidean distance \uv\ between them. For t > 1, a geometric graph G is a t- spanner, if for every pair 
of vertices u, v in V, the length of the shortest path ttg{u, v) between them in G is at most t times 
|ttu|, i.e., Vn, v G V, |vrG(n, u)| < t\uv\. A complete geometric graph on a set of points is a 1-spanner. 
Where there is no necessity to specify t, we use the term geometric spanner. A geometric spanner 
G is plane if no two edges in G cross. In this paper we only consider plane geometric spanners. A 
geometric spanner of degree at most k is referred to as a degree k geometric spanner. 

Given a geometric spanner G = {V,E), the vertex dilation or stretch factor of u,v G V, 
denoted 6g{u,v), is defined as 6g{u,v) = \7rG{u,v)\/\uv\. When G is clear from the context, 
we simply write 6{u,v). The vertex dilation or stretch factor of G, denoted 6{G), is defined as 
6{G) = sup„ v). The terms graph theoretic dilation and spanning ratio are also used in the 

literature. Refer to [22l [29l |35] for such definitions. 

Given a point set P, let the dilation of P, denoted by 5q{P), be the minimum stretch factor of 
a plane geometric graph (equivalently, triangulation) on vertex set P; see [M]. Similarly, let the 
degree k dilation of P, denoted by 6o{P,k), be the minimum stretch factor of a plane geometric 

*A preliminary version in: Proceedings of the 2nd International Conference on Algorithms and Discrete Applied 
Mathematics, Thiruvananthapuram, India, Feb. 2016, vol. 9602 of LNCS. 
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graph of degree at most k on vertex set P. Clearly, 6o{P, k) > 5o{P) holds for any k. Furthermore, 
6o{P,j) > 6o{P,k) holds for any j < k. (Note that the term dilation has been also used with 
different meanings in the literature, see for instance HUMj.) 

In the last few decades, great progress has been made in the field of geometric spanners; for 
an overview refer to |25l l35] . Common goals include constructions of low stretch factor geometric 
spanners that have few edges, bounded degree and so on. A survey of open problems in this area 
along with existing results can be found in m- Geometric spanners find their applications in the 
areas of robotics, computer networks, distributed systems and many others. Refer to miaiaiia 
for various algorithmic results. 

The existence of plane t-spanners for some constant t > 1 (with no restriction on degree) was 
first investigated by Chew [l5] in the 80s. He showed that it is always possible to construct a 
plane 2-spanner with 0{n) edges on a set of n points; he also observed that every plane geometric 
graph embedded on the 4 points placed at the vertices of a square has stretch factor at least \/2. 
This was the best lower bound on the worst-case dilation of plane spanners for almost 20 years 
until it was shown by Mulzer [3l] using a computer program that every triangulation of a regular 
21-gon has stretch factor at least (2sin^ -|-sin |y -|-sin |y )/ sin ^ = 1.4161... Henceforth, it 
was posed as an open problem by Bose and Smid m Open Problem 1] (as well as by Kanj in 
his survey [26l Open Problem 5]): ^^What is the best lower bound on the spanning ratio of plane 
geometric graphs? Specifically, is there a t > -v/2.005367532 ss 1.41611... and a point set P, such 
that every triangulation of P has spanning ratio at least t?”. We give a positive answer to the 
second question by showing that a set S of 23 points placed at the vertices of a regular 23-gon, has 
dilation So{S) > (2 sin || -|- sin ||)/ sin = 1.4308 ... 

The problem can be traced back to a survey written by Eppstein |24l Open Problem 9]; “ What 
is the worst case dilation of the minimum dilation triangulation?”. The point set S also provides 
a partial answer for this question. From the other direction, the current best upper bound of 1.998 
was proved by Xia |37! using Delaunay triangulations. Note that this bound is only slightly better 
than the bound of 2 obtained by Chew [15] in the 1980s. For previous results on the upper bound 
refer to [iniiisiiiiEHi. 

The design of low degree plane spanners is of great interest to geometers. Bose et al. |9| were the 
first to show that there always exists a plane t-spanner of degree at most 27 on any set of points in 
the Euclidean plane where t ~ 10.02. The result was subsequently improved in 13 igiTl Ea EZl I33| 
in terms of degree. Recently, Kanj et al. [28] showed that t = 20 can be achieved with degree 4. 
However, the question whether the degree can be reduced to 3 remains open at the time of this 
writing. If one does not insist on having a plane spanner. Das et al. im showed that degree 3 is 
achievable. While numerous papers have focused on upper bounds on the dilation of bounded degree 
plane spanners, not much is known about lower bounds. In this paper, we explore this direction 
and provide new lower bounds for unrestricted degrees and when degrees 3 and 4 are imposed. 

A greedy triangulation of a finite point set P is constructed in the following way: starting with an 
empty set of edges E, repeatedly add edges to E in non-decreasing order of length as long as edges in 
E are noncrossing. Bose et al. [TO] have showed that the greedy triangulation is a t-spanner, where 
t = 8(7r — afi/slv?{a/A)) k, 11739.1 and a = 7r/6. Here we obtain a worst-case lower bound of 
2.0268; in light of computational experiments we carried out, we believe that the aforementioned 
upper bound is very far from the truth. 

Related work. If Sn is the set of n vertices of a regular n-gon, Mulzer |34j showed that 
1.3836... = \J2 — \/3 -b \/3/2 < bfiSn) < 0.47l7r/sin0.47l7r = 1.4858..., for every n > 74; the 
upper bound holds for every n > 3. Amarnadh and Mitra |3| have shown that in the case of a cyclic 
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polygon (a polygon whose vertices are co-circular), the stretch factor of any fan triangulation (i.e., 
with a vertex of degree n — 1) is < 1.4846. 

As mentioned earlier, low degree plane spanners for general point sets have been studied in mm 
El na EH [33]. The construction of low degree plane spanners for the infinite square and hexagonal 
lattices has been recently investigated in |21j . 

Bose et al. [S] presented a finite convex point set for which there is a Delaunay triangulation 
whose stretch factor is at least 1.581 > '7r/2, thereby disproving a widely believed tt/2 upper bound 
conjectured by Chew m- They also showed that this lower bound can be slightly raised to 1.5846 
if the point set need not be convex. This lower bound for non-convex point sets has been further 
improved to 1.5932 by Xia and Zhang [38] . 

Klein et al. m proved the following interesting structural result. Let 5* be a finite set of points 
in the plane. Then either S' is a subset of one of the well-known sets of points whose triangulation 
is unique and has dilation 1, or there exists a number A(S) > 1 such that each finite plane graph 
containing S among its vertices has dilation at least A(S). 

Cheong et al. [14] showed that for every n > 5, there are sets of n points in the plane that do 
not have a minimum-dilation spanning tree without edge crossings and that 5 is minimal with this 
property. They also showed that given a set 5 of n points with integer coordinates in the plane and 
a rational dilation f > 1, it is NP-hard to decide whether a spanning tree of S with dilation at most 
t exists, regardless if edge crossings are allowed or not. 

Knauer and Mulzer [31] showed that for each edge e of a minimum dilation triangulation of a 
point set, at least one of the two half-disks of diameter about 0.2|e| on each side of e and centered 
at the midpoint of e must be empty of point^ 

When the stretch factor (or dilation) is measured over all pairs of points on edges or vertices of 
a plane graph G (rather than only over pairs of vertices) one arrives at the concept of geometric 
dilation of G] see [201(22] . 

Our results. (I) Let 5 be a set of 23 points placed at the vertices of a regular 23-gon. Then, 
(5o(5') = (2 sin || -|- sin ||)/sin = 1.4308 ... (Theorem Section]^. This improves the previous 
bound of (2sin^ -|-sin -|-sin |y)/ sin ^ = 1.4161..., due to Mulzer [M], on the worst case 
dilation of plane spanners. 

(II) (a) For every n > 13, there exists a set S' of n points such that <5o(S, 3) > 1 -|- \/3 = 
2.7321... (Theor em Secti on . (b) For every n > 6, there exists a set S of n points such that 

do{S, 4) > 1 -b (5 — y/h)/2 = 2.1755 ... (Theorem Section . The previous best lower bound 
of (2 sin ^ -b sin |y -b sin ^y)/ sin ^ = 1.4161..., due to Mulzer [31] holds for any degree. Here we 
sharpen it for degrees 3 and 4. 

(III) For every n > 6, there exists a set S of n points such that the stretch factor of the greedy 
triangulation of S is at least 2.0268. 

Notations and assumptions. Let P be a planar point set and G = (K, E) be a plane geometric 
graph on vertex set P. For p,q G P, pq denotes the connecting segment and \pq\ denotes its 
Euclidean length. The degree of a vertex (point) p G P is denoted by deg(p). For a specific point 
set P = {pi ,... ,pn}, we denote a path in G consisting of vertices in the order pi,pj,pk, ■ ■ ■ using 
p{i, j, k,...) and by \p{i, j, k,.. .)\ its total Euclidean length. The graphs we construct have the 
property that no edge contains a point of P in its interior. 

^ Their result inaccurately states that the entire disk of that diameter is an exclusion region. 
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2 A new lower bound on the dilation of plane spanners 


In this section, we show that the set 5 = {sq, • • •, S22} of 23 points placed at the vertices of a regular 
23-gon has dilation 5o{S) > (2 sin || + sin ||)/sin ^ = 1.4308 ... (see Fig. [^. Assume that the 
points lie on a circle of unit radius. We first present a theoretical proof showing that So{S) > (sin + 
sin || +sin ||)/ sin = 1.4237 ...; we then raise the bound to (lo('S') > (2 sin || + sin ||)/ sin = 
1.4308... using a computer program. The result obtained by the program is tight as there exists 
a triangulation of S (see Fig. [^(right)) with stretch factor exactly (2sin|| + sin||)/sin^ = 
1.4308... 


Define the convex hull length of a chord SiSj G S as = min(|i—j|, 23 —|i—j|). Observe that 

1 < hihj) — 11- Since triangulations are maximal planar graphs, we only consider triangulations 
of S while computing So{S); in particular, every edge of the convex hull of S is present. Note that 
there are C 21 = 24,466,267,020 triangulations of S. Here Cn = Catalan number 

and there are Cn ways to triangulate a convex polygon with n + 2 vertices. 

If Si,Sj G S, then |siSj| = 2sin ^^ 23 ^ ■ Consider a shortest path connecting Si,Sj G S consisting 
of k edges with convex hull lengths ni,...,nk‘, its length is \p{i,..., j)\ = 2^^^]^sin^^. Let 
A = and 


We will use A 


g{X,ni ,... ,nfc) = 




Si S i 


sin 


Att 


^zOj\ .... 23 

11 in all subsequent proofs of this section and therefore we set 


( 1 ) 


/(ni,...,nfc) := 5 (ll,ni,... ,nfc). 


( 2 ) 


Various values of /, as given by Q and Q, will be repeatedly used in lower-bounding the 
stretch factor of point pairs in specific configurations, i.e., when some edges are assumed to be 
present. Observe that / is a symmetric function that can be easily computed (tabulated) at each 
tuple ni,... ,nk‘, see Table 



Figure 1: Left: The set S of 23 points placed at the vertices of a regular 23-gon. Right: A triangulation of 
S with stretch factor (2 sin || -|- sin ||)/sin = 1.4308..., which is achieved by the detours for the pairs 
sioj S 21 and sq, sig- The shortest paths connecting the pairs are shown in blue and red, respectively. 


Given a chord SQSi, let lower(soSi) = {sj+i,..., S 22 } and upper(soSj) = {si,... Si_i}. The range 
of possible convex hull lengths of the longest chord in a triangulation of S is given by the following. 

Lemma 1. If i is the convex hull length of the longest chord in a triangulation of S, then i G 
{8,9,10,11}. 


4 









/(4,7) 

1.3396... 

❖ /(2,2,8) 

1.4308... 

/(5,6) 

1.3651... 

❖ /(3,3,5) 

1.4312... 

❖ /(5,7) 

1.4514... 

❖ /(3,4,4) 

1.4409... 

❖ /(6,6) 

1.4650... 

❖ /(1,4, 7) 

1.4761... 

/(2,3,6) 

1.4023... 

❖ /(2,3,7) 

1.4886... 

/(I,5,5) 

1.4061... 

❖ /(3,3,6) 

1.5312... 

❖ /(2,4,5) 

1.4237... 

❖ /(1,1,4, 5) 

1.4263... 

❖ /(I, 3, 8) 

1.4257... 

/(I, 2, 3, 5) 

1.4388... 


Table 1: Relevant values of f{ni,... ,nk) as required by the proofs in this section. Values used explicitly in 
the proofs are marked using ❖. 


Proof. We clearly have i > 2. Since S is symmetric, we can assume that sqSi is the longest chord. 
Since < 11 for any 0 < i,j < 22, we have I < 11. Suppose for contradiction that 2 < £ < 7. 

Then is an edge of some triangle Asos^Sm, where £ + 1 < m < 22. In particular, 

< £ < 7. (3) 

If m < 11, then = min(m, 23 — m) = m > £ + 1, a contradiction to Ps maximality. 

Assume now that m > 12; then /r(0, m) = 23 — m < i, since i is the length of a longest chord. It 
follows that m>23 — l'>23 — 7 = 16. Ifm — £<11, then /i(£, m) = m — £>16 — 7 = 9, a 
contradiction to Q. Ifm —£ > 12, then /r(£, m) = 23 —(m —£) = 23 —m + £>£ + l, a contradiction 
to £’s maximality. Consequently, we have 8 < £ < 11, as required. □ 

Proof outline. For every £ G {8,9,10,11}, if the longest chord in a triangulation T has length 
£, we show that 6{T) > /(2,4, 5) = 1.4237... Assuming that is a longest chord, we consider 
the triangle with base s-iid third vertex in upper(sos^) or lower( sqs^), depending on £. For 
each such triangle, we show that if the edges of the triangle along with the convex hull edges of 
S are present, then in any resulting triangulation there is a pair whose stretch factor is at least 
/(2,4, 5) = 1.4237... Essentially, the long chords act as obstacles which contribute to long detours 
for some point pairs. In four subsequent lemmas, we consider the convex hull lengths 8,9,10,11 
(from Lemma successively. 

In some arguments, we consider a primary pair Si,Sj, and possible shortest paths between the 
two vertices. We show that if certain intermediate vertices are present in 7r{si,Sj), then 6{si,Sj) > 
/(2,4, 5). Otherwise if certain long edges are present in 7r{si,Sj), then 6{su,Sy) > /(2,4, 5), where 
Su, Sy is a secondary pair. In the figures, wherever required, we use circles and squares to mark the 
primary and secondary pairs, respectively (see for instance Fig. [^. In some of the cases, a primary 
pair suffices in the argument, i.e., no secondary pair is needed. 

Lemma 2. //£ = 8, then 6{T) > /(2,4,5) = 1.4237... 

Proof. Refer to Fig. Let sosg be the longest chord. The triangle with base sqSs and third vertex 
in lower(soS8) has two other sides of convex hull lengths 7 and 8. It thus suffices to consider the 
triangle AsoSs^ie only and assume that the edges sosg, sgsig and so'Sie are present. 

In this proof, the primary pair is sio,S 2 i and the secondary pair is S 3 ,si 4 . Now, consider the 
pair sio,S 2 i- Note that either sq G vr(sio,S 2 i) or sie G vr(sio,S 2 i). In the former case, ( 5 (sio,S 2 i) > 
1/9(10,8, 0, 21 )|/|sioS 2 i| > /(2,8, 2) = 1.4308 ... We may thus assume that sie G 7r(sio, S 2 i)- 

Similarly, for the pair S 3 ,si 4 either so G vr(s 3 ,si 4 ) or sg G 7 r(s 3 ,si 4 ). If sq G vr(s 3 ,si 4 ), then 
< 5 ('S 3 , S 14 ) > 1/9(3,0,16, 14)|/|s3Si 4| > /(3, 7, 2) = 1.4886 ... Thus, assume that sg G '/r(s 3 , S 14 ). 
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Figure 2 : Illustrating Lemma Left: si2 G 7 r(sio,S 2 i), primary pair: sio,S2i. Right: SioSis G 7 r(sio,S 2 i), 
primary pair: Sio,S 2 ij secondary pair: 33 , 514 . 


If at least one of S 12 , S 13 , or 544 is in 7r(sio, 521 ), then 

mind^lO, 12,16, 21)|, |p(10,13,16, 21)|, |p(10,14,16, 21)|) 


<5(si0) S 21 ) > 


|'SlO'S 2 l| 


> min(/(2, 4, 5), /(3,3, 5), /(4, 2, 5)) = /(2, 4, 5) = 1.4237. 

Otherwise, one of sio^is, siosie) sii^is, or suSie must be in 7r(sio, 521 ), and 

min(|p(3,8,10,14)|, \p{3, 8,11,14)|), \p{3, 8,15,14)|) 


<^(■53,514) > 


|'S3'Sl4| 


> min(/(5, 2,4), /(5, 3, 3), /(I, 5, 7)) = /(2,4,5) = 1.4237... 
Lemma 3. If £ = 9, then S{T) > /(2,4,5) = 1.4237... 


□ 


Proof. Let sosg be the longest chord and consider the triangle with base SQSg and the third vertex 
in lower(sosg). There are three possible cases depending on the convex hull lengths of other two 
sides of the triangle: {7, 7}, {8, 6} or {9, 5}. We consider them successively. 




Figure 3 : Illustrating Case A from Lemma|^ Left: S5 G 74(33,314). Right: 3333 G 74(33,314). 
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Case A; The convex hull lengths of the other two sides are {7, 7}. Let AsoSgSie be the required 
triangle; refer to Fig. In this case, the primary pair is S 3 ,si 4 and the secondary pair is sej^is- 
Either sq G vr(s 3 , su) or sg G 7 r(s 3 , su). If sg G 7 r(s 3 , su), then 5 ( 53 , su) > |p(3, 0,16,14)|/|s3Si4| > 
/(3, 7, 2) = 1.4886 ... Thus, we assume that sg G vr(s 3 , S 14 ). 

Similarly, for the pair S 6 ,si 8 , either sg G 7r(s6,si8) or sg G 7r(s6,si8)- If sg G vr(s 6 ,si 8 ) then 
6 {sq, Sis) > |p(6,9,16,18)|/|s6Si8| > /(3, 7, 2) = 1.4886 ... Thus, assume that sg G vr(s 6 , sis). 

Now, if at least one of S 5 , sq, or sj is in 7 r(s 3 , S 14 ), then 


<I('S3,Sl4) > 


min(|p(3, 5,9,14)|, |p(3, 6 , 9,14)|, \p{3, 7,9,14)|) 


|'S3Si4| 


> min(/(2,4, 5), /(3, 3, 5), /(4, 2, 5)) = /(2,4,5) = 1.4237... 


Otherwise, one of 5388 , 5359 , 5458 , or 5459 must be in 7 r( 53 , 5 i 4 ), and 


'I('S 6 ,'S 18 ) > 


min(|/9(6,3,0,18)|, |p(6,4, 0,18)|) 

I'SgSisI 


> min(/(3,3,5),/(2,4,5)) 


/(2,4,5) = 1.4237... 


Case B: The convex hull lengths of the other two sides are {8, 6}. Let A5g595i7 be the required 
triangle; refer to Fig. [^(left). As in Case A, the primary pair is 53,514 and the secondary pair is 
56 , 518 . Consider the pair 53 , 514 . Either 547 G 77 ( 53 , 514 ) or 59 G 77 ( 53 , 514 ). If 547 G 77 ( 53 , 544 ), then 
5 ( 53 , 544 ) > 1/9(3,0,17, 14)|/|535i 4| > /(3, 6, 3) = 1.5312 ... So we assume that sg G 77 ( 53 , su)- 




Figure 4: Illustrating Case B (left) and Case C (right) from Lemma 


Similarly, for the pair 56 ,si 8 , either 59 G 7 r(s 6 , 5 i 8 ) or 5 g G 77(55, sis). H ^9 G 77(55, sis) then 
5 ( 56 , 518 ) > 1/9(6,9,17, 18 )|/|s 65 i 8 | > /(3, 8,1) = 1.4257... Thus, we assume that 5 g G 77(55,548). 
Now, it can be checked that by the same analysis as in Case A, the same lower bound of /(2,4, 5) 
holds. 

Case C: The convex hull lengths of the other two sides of the triangle are {9, 5}. Let A 5 gS 95 i 8 
be the required triangle; refer to Fig. (right). Then, 


6( 54 , 545 ) ^ 


min(|p(4,0,18,16)|, |p(4, 9,16)|) 

|'S4Sl6| 


> min(/(4, 5, 2), /(5, 7)) = /(2,4, 5) = 1.4237... □ 


Lemma 4. If i = 10, then 6{T) > /(2, 4, 5) = 1.4237... 
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Figure 5: Illustrating Case A from Lemma 


Proof. Let sqSio be the longest chord. The possible convex hull lengths of the other two sides of 
the triangle with base sqSio and the third vertex in lower(sosio) are {10, 3},{9,4},{8, 5},{7, 6}. We 
consider these cases successively. 

Case A: The convex hull lengths of the other two sides of the triangle are {10,3}, {9,4} or 
{8,5}. Let AsoSiO'S 20 ) AsqSiO'Siq, Asosio^is be the required triangles, respectively; refer to Fig. 
and Fig. [^(left). Then, 


(i(s4, srg) ^ 


min(|p(4, 0, 20,16)|, |p(4,0,19,16)|, \p{4, 0,18,16)|, |p(4,10,16)|) 


'S4'Sl6 


> min(/(4, 3, 4), /(4,4,3), /(4, 5, 2), /(6, 6)) = /(2,4, 5) = 1.4237.. 




Figure 6: Illustrating Case A (left) and Case B (right) from Lemma 


Case B: The convex hull lengths of the other two sides are {7,6}. Let AsqSioSi? be the 
required triangle; refer to Fig. [^(right). In this case, the primary pair is S 3 ,si 4 and the secondary 
pair is sei'Sis- Either sq C 7J'(s3,si 4) or sio G 7r(s3,si4). If sq G vr(s 3 ,si 4 ), then (5(s3,si4) > 
1/9(3, 0, 17, 14)|/|s3Si 4| > /(3, 6, 3) = 1.5312 ... Thus, we assume that sio G vr(s 3 , S 14 ). 

Similarly, for the pair S6,si8, either sio G 7r(s6,si8) or sq G 7r(s6,si8). If sio G 7r(s6,si8), then 
<^(•56, 'S 18 ) > 1/9(6,10,17, 18)|/|s6Si8| > /(4, 7,1) = 1.4761... Thus, assume that sq G 7r(s6, 'Si 8 ). 
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Now, if at least one of S 5 , sg, s?, or ss is in 7 r(s 3 , su), then 


<5(53, sm ) > 


min(|p(3, 5,10,14)|, |/j(3, 6,10,14)|, 1^3, 7,10,14)|, |p(3, 8,10,14)|) 


|'S3'Sl4| 


> niin(/(2, 5,4), /(3,4,4), /(4,3,4), /(5, 2,4)) = /(2,4, 5) = 1.4237 ... 
Otherwise, one of 5359 , 53510 , 5459 , or S 4510 must be in 71 ( 53 , 514 ), and 

min(|p(6, 3, 0,18)|, \p{6, 4,0,18)|) 


^(•56,518) > 


|■S65l8| 


> min(/(3, 3, 5), /(2,4,5)) = /(2,4, 5) = 1.4237... 


□ 


Lemma 5. If £ = 11, then 6 {T) > /(2,4, 5) = 1.4237... 

Proof. Let 5o5ii be the longest chord. Since the size of upper(so5ii) is smaller than the size of 
lower(so5ii), we consider upper(so5ii) is our analysis. The possible convex hull lengths of the 
other two sides of the triangle with base 5o5ii and the third vertex in upper(so5ii) are {1,10}, 
{2,9}, {3,8}, {4,7}, {5,6}. We consider the following cases successively. 




Figure 7: Illustrating Case A (left) and Case B (right) from Lemma 


Case A: The convex hull lengths of the other two sides are {1,10}. Let AsqSiSh be the 
required triangle; refer to Fig. 0(1 eft). In this case, the primary pair is 56,5i8 and the secondary 
pair is 54 , 545 . Consider the pair sg, 5i8. Either sn G 7r(s6, 5i8) or sq G 7 r(s 6 , 5i8). If 5 ii G 7r(s6, 5i8), 
then <5(56, 5i8) > |/o(6,11,18)|/|s65i8| > /(5, 7) = 1.4514... Hence, we assume that 5o G 7r(s6, sis). 

If at least one of 52 , 53 , 54 , or 35 is in 7r(s6, 513 ), then 


( 5 ( 56 , 518 ) > 


min(|p(6, 2,1,0,18)|, |p(6, 3,1,0,18)|, |p(6,4,1, 0,18)|, |p(6, 5,1,0,18)|) 


|■S65l8| 


= min(/(4,1,1,5), /(3, 2,1, 5), /(2, 3,1, 5), /(1,4,1, 5)) = /(1,1,4, 5) = 1.4263 ... 
Otherwise, 5i56 is in 7r(s6,5i8), and then 


( 5 ( 54 , 515 ) > 


min(|p(4,6,Il, 15)|, |p(4,1,0,15)|) 


54515 


> min(/(2, 5,4), /(3,1, 8)) = /(2,4,5) = 1.4237... 
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Figure 8: Illustrating Case B (left) and Case C (right) from Lemma 


Case B: The convex hull lengths of the other two sides are {2, 9}, {3, 8 } or {4, 7}. Let AsoS 2 Sii, 
AsqSsSu, AS 0 S 4 S 11 be the required triangles, respectively. Refer to Fig. [^(right) and Fig. [^(left) 
for illustrations. 

As in the Case A, we may assume that sq G vr(s 6 , sis). Then, 


<^(■56, Sis) > 


min(|p( 6 , 2,0,18)|, |p( 6 , 3, 0,18)|, \p{6, 4,0,18)|) 


Isesisl 


> min(/(4, 2, 5), /(3, 3, 5), /(2,4, 5)) = /(2,4,5) = 1.4237... 


Case C: The convex hull lengths of the other two sides are {5, 6 }. Let AsqSsSu be the required 
triangle; refer to Fig. |^(right). In this case, the primary pair is S 6 ,si 8 and the secondary pair is 
sg, S 20 - As in the Case A, we assume that sq G sig). 

Now, if at least one of sig, S 2 O) S 21 , or S 22 is in ' 7 r(s 6 , sig), then 


<^(S 6 ,S 18 ) > 


min(|p( 6 , 5,0,19,18)|, \p{6, 5,0, 20,18)|, |p( 6 , 5,0, 21,18)|, |p( 6 , 5, 0, 22,18)|) 


|S6Sl8| 


= min(/(l, 5,4,1), /(I, 5, 3, 2), /(I, 5, 2, 3), /(I, 5,1,4)) = /(1,1,4, 5) = 1.4263 ... 


Otherwise, sqSis is in ' 7 r(s 6 , sig), and then 

min(|p(8,5,0,20)|,|p(8,ll,18,20)|) 


5(s8,S20) > 


|S8S20| 


> min(/(3, 5, 3), /(3, 7, 2)) = /(3, 3, 5) = 1.4312 . 


□ 


Putting these facts together yields the main result of this section: 

Theorem 1. Let S be a set of 23 points placed at the vertices of a regular 23-gon. Then 

. llvr 


^o{S) = f{2, 2 , 8 )= ( 2 sin ^ + sin 


23 


23 J 


sm-= 1.4308... 

23 


Proof. By Lemmas 2][5, we conclude that 6o{S) > f{2, 4, 5) = (sin || + sin || + sin ||)/sin ^ = 
1.4237... On the other hand, the triangulation of S in Fig. [^(right) has stretch factor /(2, 2, 8) = 
1.4308 ... and thus /(2,4, 5) = 1.4237 ... < 5o{S) < /(2, 2,8) = 1.4308 ... 
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A parallel C++ progranij^ that generates all triangulations of S based on a low memory algo¬ 
rithm by Parvez et al. |36[ Section 4] shows that each of the C 21 triangulations has stretch factor 
at least /(2,2,8). We thereby obtain the following final result: So{S) = /(2,2,8) = (2sin|| -|- 
sin ||) / sin = 1.4308 ... □ 

Remarks. Using the program we have also checked that the next largest stretch factor among 
all triangulations is /(3,3, 5) = 1.4312..., and further that there is no triangulation of S that has 
stretch-factor < 1.4312 other than /(2,2,8). Thus, the result in Theorem is not affected by 
floating-point precision errors. 



S? (S6 


S5 



Figure 9: Triangulations of 825 and S '26 with stretch factors < 1.4296 and < 1.4202, respectively. Worst 
stretch factor pairs are marked in circles and the corresponding shortest paths are shown in red. 


Let Sn denote the set of points placed at the vertices of a regular n-gon. Using a computer 
program, Mulzer obtained the values <5o(*S'n) for 4 < n < 21 in his thesis [34t Chapter 3]. Using our 
C++ program, we confirmed the previous values and extended the range up to n = 24: 5 q{S 22 ) = 
1.4047..., So{S 24 :) = 1.4013... and somewhat surprisingly, 60 ( 323 ) = 1.4308... By upper bound 
constructions, it follows that (5o(5'25) < 1.4296 and 6 o{S 2 o) < 1.4202; see Fig.Observe that 6 o(Sn) 
does not exhibit a monotonic behavior; see Table 


n 

6 o{Sn) 

n 

<5o(5n) 

n 

5o(5n) 

4 

1.4142... 

12 

1.3836... 

20 

1.4142... 

5 

1.2360... 

13 

1.3912... 

21 

1.4161... 

6 

1.3660... 

14 

1.4053... 

22 

1.4047... 

7 

1.3351... 

15 

1.4089... 

23 

1.4308 ... 

8 

1.4142... 

16 

1.4092... 

24 

1.4013... 

9 

1.3472... 

17 

1.4084... 

25 

< 1.4296 

10 

1.3968... 

18 

1.3816... 

26 

< 1.4202 

11 

1.3770... 

19 

1.4098... 




Table 2: The values of So(Sn) for n = 4,..., 26. 


^Refer to the .cpp file and the Appendix within the source at arXiv: 1509.07181 
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3 Lower bounds for the degree 3 and 4 dilation 

In this section, we provide lower bounds for the worst case degree 3 and 4 dilation of point sets in 
the Euclidean plane. We begin with degree 3 dilation. We first present a set P of n = 13 points (a 
section of the hexagonal lattice with six boundary points removed) that has 6o{P, 3) > 1 + ^/3 and 
then extend P to achieve this lower bound for any n > 13. 

Theorem 2. For every n > 13, there exists a set S of n points sueh that 6 q{S,3) > 1 + \/3 = 
2.7321... The inequality is tight for the presented sets. 

Proof. Let P = {po} U Pi U P 2 be a set of 13 points as shown in Fig. [^(left) where Pi = 
{pi,P 3 ,P 5 ,P 7 ,P 9 ,Pn} and P 2 = {p 2 ,Pi,P 6 ,Ps,Pio,Pi 2 }- The points in Pi and P 2 lie on the vertices 



Figure 10: Left: the point set P = {po,pi,... ,pi 2 }; some pairwise distances are: IP 2 P 12 I = 2, \p 2 P 3 \ = 
IpiPsI = IP 1 P 12 I = 's/S. Right: a plane degree 3 geometric spanner on P with stretch factor 1 + -s/S, which is 
achieved by the detours for the point pairs {pijPs}, {P 5 ,P 7 } and {p 9 ,Pii}. 


of two regular homothetic hexagons centered at po of radius 1 and 2 respectively. Furthermore, the 
points in each of the sets {p 2 ,Pi,Po,P 7 ,P 8 }, {P 4 .,P 3 ,Po,P 9 ,Pio} and {pi 2 ,Pii,Po,P 5 ,P 6 } are collinear. 

We show that 6o{P,3) > 1 + \/3. Since no edge can contain a point in its interior, the point 
Po can have connecting edges only with the points from Pi. First, assume that the six edges in 
E = {P 1 P 3 ,P 3 P 5 ,P 5 P 7 ,P 7 P 9 ,P 9 P 11 ,P 1 P 11 } are present (see Fig. [IT](left)). 



•P8 PG- 'PS PG' 

Figure 11: Left: all edges in E are present. Right: Case B. 


We can also assume that the edge popi is present since po must be connected to at least one of 
the points in Pi. Observe that now deg(pi) = 3. In this case. 


d{pi,P2) > 


\PlP2\ 


> 1 + \/3, where i G {3,11}. 
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Now assume that an edge in E, say piPs, is missing. Then, the following three cases arise depending 
on deg(po) G {1,2,3}. 

Case A: If deg(po) = 1, then 


5(P1,P3) > 


1 /^( 1 , ^, 3)1 

IpiPsI 


> 1 + \/3 where i G (2, 5,4,11}. 


Case B: If deg(po) = 2, consider the edges PoPi,PoP 3 ', see Fig. 11 (right). If PoPi,PoP 3 are 
present 6 {po,py) > |/9(0, 3, 7)|/|poP7| = 1 + \/3 else if at least one edge in {poPi,PoP 3 } is absent then 
since pip 3 is absent, 6 {pi,p 3 ) > 1 + \/3 by the same analysis as in Case A. 

Case C; If deg(po) = 3, then if at least one of the edges PoPi,PoP 3 is absent, S{pi,p 3 ) > l + \/3 as 
shown in CASE A. Thus, assume that poPi,poP 3 are present. Now, the following two non-symmetric 
cases will arise. Either pop^ is present or poPY is present. 


P12 - 


• P2 


Pl2 . 


• P2 



Figure 12: Left: the edges PoPi,PoP 3 ,PoP 5 are present. Right: the edges PoPi,PoP 3 ,PoP 7 are present. 


If PoP5 is present (refer to Fig. [^(left)) then. 


^{po,P 9 ) > > 1 + V3, where i G {1,5}. 

\P0P9\ 


Now assume that poP 7 is present (refer to Fig. 12 (right)). Observe that if pyP 9 is absent then, 

|p(7,7,9)1 


^{P7,P9) > 


\P7P9\ 


> 1 + V3, where i G {8,10,11}. 


Thus, assume that PYP 9 is present. Similarly, assume that p^pY is present, otherwise 

1/3(7,7,5)1 


KP7,Pb) > 


\P7Pb\ 


> 1 + Vs, where i G {3, 5, 6}. 


Now, as PoPy, P 5 Py and pyPq are present, deg(p 7 ) = 3. In this case, 

VP 7 ^pV = > 1 + Vs, where i G {5,9}. 

\P7P3\ 

We have thus just shown that <50(73,3) > 1 + \/3. For n > 14, we may assume that po = (0,0), 
P 3 = (1,0), and let pi = (x + i, 0) for 1 = 13,..., n — 1, where x 3> 1 (e.g., setting x = 100 suffices); 
finally, let S = PUP', where P' = {pis,... ,Pn-i}. If 7t,n G P C S', then going from n to n via P' is 
inefficient, so as shown earlier in this proof, 6 {u, x) > 1 + VS. Thus, <5o(S, 3) > 1 + VS, as required. 
Moreover, this lower bound is tight for both P and S; see Fig. 11 (right). □ 
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Remark. If A is the infinite hexagonal lattice, it is shown in |21] that 3) = 1 + ^/3. 

We now continue with degree 4 dilation. We first exhibit a point set P of n = 6 points with 
degree 4 dilation 1 + \J and then extend it so to achieve the same lower bound for any 
larger n. Consider the 6 -element point set P = {po ,... ,^ 5 }, where pi,... ,p 5 are the vertices of a 
regular pentagon centered at po- 


Theorem 3. For every n>6, there exists a set S of n points such that 

5oiS, 4) > 1 -h = 2.1755... 

The inequality is tight for the presented sets. 


Proof. Assume that pi,... ,p 5 lie on a circle of unit radius centered at po- Since deg(po) < 
exists a point p*, 1 < i < 5 such that pop* is not present; we may assume that i = 1; see 
Observe that 


4, there 


Fig. 13 



Figure 13: A plane degree 4 geometric graph on the point set {po, ■ ■ ■ jPs} that has stretch factor exactly 
1 -I- y (5 — •\/5)/2, which is achieved by the detour between the pair poiPi- 


IpoPiI = 1 and \piP 2 \ = IpiPsI = v^l^ + 1 ^ - 2 • 1 • 1 cos(27r/5) = yj(5 - V5)/2. 


Now, 


^iPo,Pi) > 


1/3(0,bl)| 

\PoPi\ 


> 1 -b y^(5 - V5)/2 = 2.1755 ..., where i E { 2 , 5}. 


Thus, (5o(P, 4) > 1 + y^( 5 — \/5)/2. As in the proof of Theorem the aforesaid six points can 
be used to obtain the same lower bound for any n > 6. 


To see that the above lower bound is tight, consider the degree 4 geometric graph on P in Fig. 13 
whose stretch factor is exactly that, due to the detour between po,pi. □ 


4 A lower bound on the dilation of the greedy triangulation 

In this section, we present a lower bound on the worst case dilation of the greedy triangulation. 
Place four points at the vertices of a unit square U, and two other points in the exterior of U on 
the vertical line through the center of U and close to the lower and upper sides of U, as shown in 
Fig. [^(left). For any small e > 0, the points can be placed so that the resulting stretch factor is at 
least (J(po,P 3 ) > 2 — e. A modification of this idea gives a slightly better lower bound. 

Theorem 4. For every n > 6, there exists a set S of n points such that the stretch factor of the 
greedy triangulation of S is at least 2.0268. 
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Proof. Replace the unit square by a parallelogram V with two horizontal unit sides, unit height and 


angle a G (vr/4,7r/2) to be determined, as shown in Fig. 14 (right). Place four points at the vertices 
of V and two other points in the exterior of V on the vertical line through the center of the V and 
close to the lower and upper side of V. First, observe that the greedy triangulation is unique for 
this point set. Second, observe that there are two candidate detours connecting po with p^: one of 
length (slightly longer than) 1 + a and one of length (slightly longer than) 2x + b, where a is the 
length of the slanted side of P, 6 is the length of the short diagonal of V, and x is the horizontal 
distance between the upper left corner and the center of V. 


Po 



Pi 


P2 



Figure 14: Greedy triangulation of 6 points with stretch factors 2 — e (left) and 2.0268 (right). 
A straightforward calculation gives: 


1 , \/l + sin^ a — 2 sin a cos a . 1 — cot a 

a = ^-, 0 --^ and x =-. 


sm a 


sina 


T . ^ -A 1 sin^ a — 2 sin a cos a \ /tt vr 

Let fia) = mm 1 H-, 1 — cot a H-, for a G —, — 

\ sm a sm a j V 4 2 

Setting a = 1.3416 (i.e., a = 76.87°) yields 

^{Po,P3) > max /(a) >/(1.3416) = 2.0268 ..., 

Q:G('7r/4,7r/2) 

as required. As in the proofs of Theorems and the lower bound can be extended for every n > 6 
in a straightforward way. □ 


5 Concluding remarks 

In Section we have shown that any plane spanning graph of the vertices of a regular 23-gon 
requires a stretch factor of (2 sin ^ + sin ||)/ sin ^ = 1.4308 ... Henceforth, the question of Bose 
and Smid m Open Problem 1] mentioned in the Introduction can be restated: 

Problem 1. Does there exist a point set S in the Euclidean plane such that 6o{S) > (2sin^ + 
sin ||) / sin ^ = 1.4308 ...? 

Next in Section it has been shown that there exist point sets that require degree 3 dilation 

1 + -\/3 = 2.7321... (Theoremj^ and degree 4 dilation 1 + \J— 'Jb)/2 = 2.1755... (Theorem j^. 
Perhaps these lower bounds can be improved. 
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Problem 2. Does there exist a point set in the Euclidean plane that has degree 3 dilation greater 
than 1 + "v/S ? Does there exist a point set in the Euclidean plane that has degree 4 dilation greater 

than 1 + y^(5 — V^)l2 ? 

Finally in Section]^ we show that the stretch factor of the greedy triangulation is at least 2.0268, 
in the worst case. Perhaps this lower bound is not far from the truth. Using a computer program we 
have generated 1000 random uniformly distributed n-element point sets in a unit square for every 
n in the range 4 < n < 250, and computed the greedy triangulations and corresponding stretch 
factors. The highest stretch factor among these was only 1.97 (as attained for a 168-element point 
set), and so this suggests the following. 


Problem 3. Is the worst case stretch factor of the greedy triangulation attained by points in convex 
position? 


Observe that the point set used in the lower bound construction in Theorem is convex, so it 
is natural to ask: given a non-convex point set S and a greedy triangulation of S having stretch 
factor A, does there always exist a convex subset S' C S such that the stretch factor of a greedy 
triangulation for S' also equals A? The point set S = {pi,... ,pq} illustrated in Fig. 15 shows that 
this is not the case. It is routine to verify that the stretch factor of the greedy triangulation of each 
convex subset S' C S is at most 1.4753... < A = 1.4772 ... 



Figure 15: Left: greedy triangulation of a set of 6 points not in convex position with stretch factor A = 

1.4772.. . attained by the pair {pi,p 4 }. Right: the largest stretch factor of the greedy triangulation of a 
convex subset is that for the subset S' = {pi,p 2 ,P 4 ,P 6 }', h is attained by the same pair {pi,P 4 } and equals 

1.4753 .. . < A. The corresponding shortest paths are drawn in red color. 
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Appendix 


Source code. The following parallel C++ code is written using OpenMP in C++11 (notice the 
pragma directives present in the code). For the set of = 23 points, the program ran for ap¬ 
proximately 2 days on a computer with quad core processor. The program was compiled with g++ 
4.9.2. Please note that older versions of g++ might have issues with OpenMP support. Following 
is a correct way of compiling the program. 


g++ program.cpp -std=c++ll -fopenmp -03 


The number of threads has been set to 4 using the variable numberOfThreads in main(). The user 
may alter the value of the variable depending on the processor. 

Following is the output from the program. 


Execution started... 

Triangulations checked: 24466267020 
Dilation: 1.4308143191 

Time taken: 162829 seconds 


#include <iostream> 
ttinclude <cmath> 

#include <omp.h> 

#include <list> 

#include <vector> 

#include <chrono> 

#define M_PI 3.14159265358979323846 

using namespace std; 
using namespace chrono; 

struct Edge { unsigned u,v; }; 

struct gc { unsigned vis_v; Edge oppositeEdge; I; 
unsigned numberOfThreads; 

vector<unsigned long long> countTriangulations; 
vector<double> minStretchFactor; 
list<unsigned> jobList; 

typedef vector<vector<double>> Matrix2D; 
vector<Matrix2D> distarrayCollection; 

struct Triangulation 

{ 

unsigned N; 
list<gc> gcList; 
list<Edge> bcList; 
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Triangulation (const unsigned numberOfPoints) 

N = numberOfPoints; 

for(unsigned vertexID = 2; vertexID < N-1; vertexID++) 

gcList.push_frent({vertexID,{vertexID-1,vertexID+1}}); 

> 

void flipgc (const unsigned gchord) 

{ 

list<gc>::iterator hold; 

for(auto it = gcList .beginO ; it != gcList. endO ; it++) 
if(it->vis_v == gchord) 

{ 

bcList.push_front(it->oppositeEdge); 

if (next (it, 1) != gcList. endO ) 

next (it,l)->oppositeEdge.V = it->oppositeEdge.v; 

if (it != gcList .beginO ) 

prev(it ,l)->oppositeEdge.u = it->oppositeEdge.u; 

hold = it; 
break; 

> 

gcList. erase (hold); 

> 

>; 

inline double distance (const unsigned pi, const unsigned p2, const unsigned N) 

{ 

unsigned absVal = max(pl,p2) - min(pl,p2); 
unsigned lambda = min(absVal,N-absVal); 
return 2*sin((lambda*M_PI)/N); 

> 


void calculateStretchFactorDfTriangulation(const Triangulation T, const unsigned 

thread_ID) 

{ 

double stretchFactor = 0; 

Matrix2D dist = distarrayCollection[thread_ID]; 

for(unsigned i = 0; i < T.N; i++) 

for(unsigned j = 0; j < T.N; j++) 

dist[i][j] = (i != j)? INFINITY : 0; 

for(unsigned i = 0; i < T.N-1; i++) 

dist[i][i+l] = dist[i+l][i] = distance(i,i+l,T.N) ; 

dist [T.N-1][0] = dist[0][T.N-1] = distance(T.N-l,0,T.N); 

for(auto it = T.gcList.beginO ; it != T. gcList. endO ; it++) 

dist[0][it->vis_v] = dist[it->vis_v] [0] = distEince(0,it->vis_v,T.N); 
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for(auto it = T.bcList.beginO ; it != T.bcList.endO ; it++) 

dist[it->u] [it->v] = dist[it->v][it->u] = distEince(it->u,it->v,T.N); 

for(unsigned k = 0; k < T.N; k++) 

for(unsigned i = 0; i < T.N; i++) 

for(unsigned j = 0; j < T.N; j++) 

if (dist[i][j] > dist[i][k] + dist[k][j]) 
dist[i][j] = dist [i] [k] + dist[k][j]; 

for(unsigned i = 0; i < T.N; i++) 

for(unsigned j = 0; j < T.N && i != j; j++) 

{ 

double tempRatio = dist[i][j] / distan.ce(i, j ,T.N) ; 
if (tempRatio > stretchFactor) 
stretchFactor = tempRatio; 

} 

if (minStretchFactor[thread_ID] > stretchFactor) 
minStretchFactor[thread_ID] = stretchFactor; 

} 

void findChildren(Triangulation &T, const unsigned gchord, const unsigned thread_ID) 

{ 

T.flipgc(gchord); 
countTriaugulations[thread_ID]++; 

calculateStretchFactorOfTriangulation(T,thread_ID); 
if (!T.gcList.empty 0) 

for(auto it = T.gcList .beginO ; it != T.gcList .endO ; it++) 
if(it->vis_v >= T.bcList.front0.u) 

{ 

Tricuigulation *childTricLngulation = new Triaiigulation(T) ; 
findChildren(*childTriangulation,it->vis_v,thread_ID); 
delete childTriangulation; 

> 

} 

void thread_job (const unsigned N, const unsigned thread_lD) 

{ 

Triangulation *childTriangulation; 
unsigned gchord; 

while (true) 

{ 

gchord = 0; 

#pragma omp critical 

{ 

if (jobList. sizeO > 0) 

{ 

gchord = jobList.front 0; 
jobList.pop_front(); 

> 
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} 


ifCgchord == 0) 
break; 

childTriangulation = new TriELngulation(N); 
findChildren(*childTriangulation,gchord,thread_ID); 
delete childTriangulation; 

> 

} 

void findAllTriangulations (const unsigned N) 

{ 

Triangulation rootTricUigulation(N); 
countTriangulations[0]++; 

calculateStretchFactorOfTriangulation(rootTricUigulation,0); 

#pragma omp parallel for num_threads(numberOfThreads) 
for(unsigned thread_id = 0; thread_id < numberOfThreads; thread_id++) 
thread_job(N,thread_id); 


} 

void initDataStructures (const unsigned N) 

minStretchFactor.assign(numberOfThreads,INFINITY); 
countTriangulations.assign(numberOfThreads,0); 

for(unsigned gc = N-2; gc >= 2; gc—) 
j obList.push_frent(gc); 

distarrayCollection.resize(numberOfThreads); 

for(unsigned thread_ID = 0; thread_ID < numberOfThreads; thread_ID++) 

{ 

Matrix2D matrix(N); 

for(unsigned pos = 0; pos < N; pos++) 
matrix[pos].resize(N); 
distarrayCollection[thread_ID] = matrix; 

> 

} 

int mainO 

{ 

unsigned N = 23; 
numberOfThreads = 4; 

double dilation = INFINITY; 

unsigned long long totalNoOfTrieuigulations = 0; 

cout << "Execution started...\n" ; 

auto start = system_clock::now(); 
initDataStructures(N); 
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findAllTriangulations(N); 

for(unsigned thread_ID = 0; thread_ID < numberOfThreads; thread_ID++) 

{ 

if ( countTriangulations[thread_ID] > 0 && minStretchFactor[thread_ID] < dilation) 
dilation = minStretchFactor[thread_lD]; 
totalNoOfTriangulations += countTriangulations[thread_lD]; 

> 

printf ("Triangulations checked; "/olld\n" .totalNoOfTriangulations) ; 
printf ("Dilation: %1.lOf\n" .dilation); 

system_clock::time_point stop = system_clock::now(); 

auto duration = duration_cast<seconds>( stop - start ).count(); 

cout « "Time taken: " << duration << " seconds" « endl; 

return EX1T_SUCCESS; 
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